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Abstract 

We consider a quantum system composed of a spatially infinitely 
extended free Bose gas with a condensate, interacting with a small 
system (quantum dot) which can trap finitely many Bosons. Due to 
spontaneous symmetry breaking in the presence of the condensate, 
the system has many equilibrium states for each fixed temperature. 
We extend the notion of Return to Equilibrium to systems possess- 
ing a multitude of equilibrium states and show in particular that a 
condensate coupled to a quantum dot has the property of Return to 
Equilibrium in a weak coupling sense: any local perturbation of an 
equilibrium state of the coupled system, evolving under the interact- 
ing dynamics, converges in the long time limit to an asymptotic state. 
The latter is, modulo an error term, an equilibrium state which de- 
pends in an explicit way on the local perturbation (an effect due to 
long-range correlations). The error term vanishes in the small cou- 
pling limit. 

We deduce the stability result from properties of structure and 
regularity of eigenvectors of the generator of the dynamics, called the 
Liouville operator. Among our technical results is a Virial Theorem 
for Liouville type operators which has new applications to systems 
with and without a condensate. 
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1 Introduction 



This paper is devoted to the study of the dynamics of a class of quantum 
systems consisting of a small part in interaction with a large heat reservoir, 
modelled by an infinitely extended ideal gas of Bosons. We further develop 
spectral methods in the framework of algebraic quantum field theory and 
apply them to the class of systems at hand, for which the already existing 
techniques cannot be applied. 

Our main physical interest is the long-time behaviour of initial states 
close to an equilibrium state of a Bose gas that is so dense (for fixed temper- 
ature) or so cold (for fixed density) that it has a Bose-Einstein condensate, 
in interaction with a small system that can store a finite number of Bosons. 
The condensate induces long-range correlations in the system and as a conse- 
quence there are many equilibrium states for a fixed value of the temperature 
T = 1/ (3. It is thus necessary to extend the notion of Return to Equilibrium 
to systems possessing a multitude of equilibrium states. In general a system 
has a special class of equilibrium states (extremal ones) whose superposi- 
tions generate all equilibrium states. It is reasonable to expect (and proven 
for the model considered here) that each extremal equilibrium state has the 
property of Return to Equilibrium. This leads to a general definition of this 
property which we introduce in Section 1.1.2. A feature of this situation is 
that starting with a local perturbation of a given superposition of extremal 
equilibrium states the system converges in the long time limit to a possi- 
bly different superposition of extremal equilibrium states (redistribution of 
phases). The asymptotic state depends thus on the initial condition. This 
limitation of the dispersive nature of a system is natural in view of the long- 
range correlations. 

One of our goals is to prove weak coupling return to equilibrium, in the 
above setting, saying that any initial condition close to an equilibrium state 
of the coupled system, evolving under the coupled dynamics, converges in the 
long time limit to an asymptotic state. The latter is again an equilibrium 
state (different from the initial one), modulo an error which disappears in 
the hmit of small couphng. We expect a stronger result to hold, namely that 
any initial condition close to the interacting system converges, in the long 
time limit, to an equilibrium state of the interacting system. This result has 
been obtained for systems without a condensate in a variety of recent papers, 
[JP1,BFS,M1,DJ,FM2]. It is surprising that none of the methods developed 
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in these references - nor elsewhere, according to our knowledge - can be 
applied to the present case. This is due to the fact that the form factor 
of the interaction, a coupling function g e L^(R^, (Pk), whose properties are 
dictated by physics, exhibits the infrared behaviour < |g'(0)| < oo. It lies in 
between the two "extreme" behaviours g{0) =0 (more precisely, g{k) ~ \k\^, 
some p > 0, as \k\ ~ 0) and |5'(0)| = oo (more precisely, g{k) ~ as 
\k\ ~ 0), which are the only ones that can be treated using approaches ex- 
isting so far. We give here a partial remedy to this situation by establishing 
a "positive commutator theory" (a first step in a Mourre theory) which is 
applicable to a wide variety of interactions, including the case where g{0) is 
a nonzero, finite constant. Our remedy is only partial in that so far, we show 
that the equilibrium state is stable, in the sense mentioned above, but we 
cannot prove return to equilibrium. The obstruction seems to be of technical 
nature (see Section 2.2.1 for a discussion of this point). 

After clarifying (defining) the notion of Return to Equilibrium in the set- 
ting of multiple KMS states our analysis consists of two main steps. The 
first one is to carry out the decomposition of an arbitrary equilibrium state 
into extremal ones and to describe the dynamics of each of them. The sec- 
ond step in our analysis, which is the main part of this paper, consists in 
analyzing the time asymptotic behaviour of each extremal equilibrium state. 
We do this by examining the spectrum of the Liouville operators generating 
the dynamics. Our approach gives an extension of the positive commutator 
method, including a new virial theorem which yields improved results in the 
analysis of related problems for systems without a condensate. 

1.1 An extended notion of Return to Equilibrium 

We review the notion of Return to Equilibrium and extend it to systems 
having multiple KMS states. Our guiding example is a reservoir of free 

Bosons where the non-uniqueness of KMS states is due to spontaneous gauge- 
group symmetry breaking in the presence of a Bose-Einstein Condensate. 

1.1.1 Free Bose gas and Bose-Einstein Condensation 

The kinematical algebra describing the Bose gas is the Weyl algebra W{V) 
over a suitably chosen test-function space of one-particle wave functions V C 
L'^{R^,d^k). W{V) is generated by Weyl operators W{f), f eV, satisfying 
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the canonical commutation relations (CCR) 

W{f)W{g) = e-^^-<^'^)iy(/ + g), (1) 

where (•, •) is the inner product induced by L^(M^, d^/c). The dynamics of 
the Bose gas is given by the Bogoliubov transformation 

W{f) ^ at{W{f)) = W{e''^f), (2) 

where 

uj{k) = \k\^, or uj{k) = \k\. (3) 

The first choice in (3) describes non-relativistic Bosons, while the second one 
describes massless relativistic ones. 

We outline first the construction of the equihbrium state of [AW] , which 
gives a good physical understanding of the emergence of a condensate, and 
then we relate this to the works of [C] and [LP]. Any state uj on the Weyl alge- 
bra W(V) is uniquely determined by its so-called generating (or expectation) 
functional E : V ^ C, given by 

u;{W{f))^Eif), (4) 

and conversely, if : D — > C is a (non linear) function satisfying certain 
compatibility conditions then it defines uniquely a state on W{'D), see e.g. 
[A, M2]. 

Let 9 A; 1-^ p{k) > be a given function (the "continuous momentum- 
density distribution"), and po > a fixed number (the "condensate density"). 
Araki and Woods [AW] obtain a state of the Bose gas by the following pro- 
cedure. Restrict the gas to a periodic box A of volume V in and put 
Vpo particles in the ground state of the one particle Hamiltonian H\ = —A 
(or H\ = \/—A), and a discrete distribution of particles in excited states. 
Then take the limit V ^ oo while keeping po fixed and letting the discrete 
distribution of excited states tend to p{k). Like this [AW] obtain the family 
of generating functionals 



<:j(/)=exp 



(l + 2(27r)V) /> 



Jr 



(V'2(2^|/(0)|) , (5) 



where Jo{^/a'^ + (3^) = J^^ £e-^("^°"^+'5'^^°^), a,(3 e R (Bessel function). 
E^^^ defines uniquely a state of the infinitely extended Bose gas according 
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to (4) . The physical interpretation is that the resuhing state describes a free 
Bose gas where a sea of particles, all being in the same state (corresponding 
to the ground state of the finite- volume Hamiltonian) , form a condensate 
with density po, which is immersed in a gas of particles where p{k) particles 
per unit volume have momentum in the infinitesimal volume d^k around 
A; G M'^. Since the Hamiltonian in the finite box is taken with periodic 
boundary conditions the condensate is homogeneous in space (the ground 
state wave function is a constant in position space). The resulting state is 
an equilibrium state (satisfies the KMS condition) if the momentum density 
distribution is given by 

pik) = i^^r'^^^j^,^ (6) 

corresponding to Planck's law of black body radiation. 

According to the principles of quantum statistical mechanics, the equilib- 
rium state of the infinite system is obtained by taking the thermodynamic 
limit of local canonical or grand-canonical Gibbs states. This has been done 
in [C] (canonical) and [LP] (grand-canonical). We first outline the result of 
[LP]. The density matrix (acting on Fock space) for the local system is 

where z = e^^ is the fugacity and A'a is the number operator. For a fixed 
inverse temperature < /3 < oo define the critical density by 

PcHt(/3) = (27r)-|-^, (8) 

and denote by p the total (mean) density of the gas. Then 

EtU) tr {ai^WU)) Ep-,{fl (9) 

where the hmit lAI — > oo is taken with the constraint p = -t-ttt^^—t^, deter- 
mining the value of z{K). The limiting generating functional (9) is 

E,-Jf) = I e-m\<^'^^^), -p < pUP) (10) 
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where, with pQ = p- pcrit(/3) > 0, 



2 

(11) 

and p = p{k) is given in (6). For subcritical density, the number 2:00 G [0, 1] 
is determined by the equation 



(2^)-' / ZK^d'k. (12) 



00 



In the supercritical case we have z^o — 1 which corresponds to a vanishing 
chemical potential, p^ = 0. 

The thermodynamic limit of the canonical local Gibbs state is treated in 
[C], the density matrix is 



^ g -I3Ha 

^^'^ ^ tre-P^-Py ^^^^ 



and Ppv is the projection onto the subspace of Fock space with pV parti- 
cles (if pV is not an integer take a convex combination of canonical states 
with integer values piV and P2V extrapolating pV). The limiting generating 
functional is given by 

. .^ [ e-ill/ll'e~H^'^^^>, p < p,Ht(/5) (14) 

It coincides with the grand-canonical generating functional in the subcritical 
case, and with the Araki- Woods generating functional with p given by (6) 
and 

PO^ P- Pcrit 

in the supercritical case. To complete the exposition of this triangle of gen- 
erating functionals we mention that (see [C]) the grand-canonical and the 
canonical generating functionals are linked in the supercritical case, po > 0, 
by the Laplace transform 

E^)- / K{r;p)El(f)dr, (15) 
Jo 
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where the Kac density K{r; p) is 

,16, 

This means that the grand-canonical equihbrium state with supercritical 
mean density p is a superposition of canonical equilibrium states with su- 
percritical densities r, weighted with the Kac density K{r,p). 

To ease the notation we simply write uup for the equilibrium state corre- 
sponding to (15) (imagining a supercritical mean density p to be fixed). Space 
translations are given by Txf{y) — f{y — x), for a; E M.^ {Txf{k) — e'^'^^f{k) 
in Fourier space) and since E^°^{Txf) = E^°^{f) for all a; e (see (11)) the 
state ujj3 is space translation invariant. However, due to the presence of the 
condensate, the system has long range correlations: 

hm u;fs{W{f)W{rxg)) 

|a;|— >oo 

= ujp{W{f )) co^iWig)) exp [Sn'poRe (7(0)(7(0))] , (17) 

where po = p — pcrit is the condensate density. This means that cj/j is not a 
factorial state, i.e. the von Neumann algebra of observables represented in 
the Hilbert space associated to {W{V),ujp) is not a factor (see for instance 
[Ha], Theorem 3.2.2). In Section 2 we decompose cup into a superposition of 
extremal factorial states 

u;p{A) = / dpf,-,{04{A), (18) 
where the probability measure 

df) 

dpp-,iO -.^ K{r,p)dr- (19) 

is supported on = (r, 9) G [pcrit, oo) x S^} C M^. Each c<j| is a /j-KMS state 
w.r.t. the dynamics (2), having the cluster property (compare with (17)!) 

hm uj${W(f)W(Txg)) = 4{W(f)) u;${W(g)), (20) 

which is also called the property of strong mixing w.r.t. space translations. 
Clearly, the gauge transformations 7s(H^(/)) = W{e^^f), s e M, commute 
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with the dynamics given in (2), ato^is = Js°(^t for all s, t G M, and the state uj/3 
is invariant under 7^. However, the equilibrium states cj| are not invariant 
under 7^; hence there is a family of equilibrium states possessing "less 
symmetry than the dynamics" , a property of the system called spontaneous 
symmetry breaking. 

We can take any probability measure /i, supported inside [pcrit, 00) x S^, 
and define the /3-KMS state 



This shows that there is, for each (3 fixed, a multitude of equilibrium states 
(each determined by a /i), and we must examine what this means for the 
notion of Return to Equilibrium. 

Remark. The integration over r > Pcrit in (18) comes from the use 
of the grand-canonical ensemble, see (15). In the canonical case the value 
of r is fixed and the integration is only over 9 e S^. Each state lo^ is 
extremal invariant for space translations {uj^ is r^^-invariant and cannot 
be decomposed into a convex combination of two other Tx-invariant states; 
this follows from the cluster property (20)). As a consequence the average 
lim|A|^oo jx] I A '^^•^ ^?(''"^(^)) exists and is a multiple of the identity (^4 a local 
observable). In particluar, for each ^ fixed, the density p and the phase 9 are 
observables in 7r^(2t)" which reduce to numbers, (p is obtained by choosing 
Tx{A) = n{x) = a*{x)a{x) representing the number of particles at the point 
X, e*^ is obtained for Tx{A) = a{x), see [H].) The density and the phase are 
observables which commute with all other observables (they belong to the 
center of the observable algebra associated with cujj) and are thus regarded 
as classical quantities. For instance, a maximal lack of knowledge of the 
value of 9 is expressed by taking a superposition of uj^ where 9 is uniformly 
distributed. In the state a;| the values of density and phase are fixed. 

1.1.2 Return to Equilibrium 

Consider first the equilibrium state Ujs describing a reservoir of free Bosons 
at inverse temperature /? with a fixed subcritical density p < Pcrit (no con- 
densate). This system has the property of Return to Equilibrium: for any 




(21) 
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observables A,B we have 

lim ujp{B*at{A)B) = ujf3{B*B)up{A). (22) 

The functional A i-^ ujf^{B*AB) defines a vector state (if ujg{B*B) = 1) 
which is a local perturbation of LUp. Equation (22) extends to the closure 
of all convex combinations of such vector states (called the set of uj p-normel 
states). Property (22) has been shown to hold also for reservoirs (without 
condensates) interacting with a small system; then in (22) is the equi- 
librium state of the interacting system and at is the coupled dynamics, see 
[JP1,BFS,DJ,M,FM2]. For different (scattering) approaches to similar prob- 
lems we refer to [R,LV], and to [HL,QV] for stochastic methods. 

Consider next the Bose gas is in a state cup with a condensate. As we 
have seen above there are many (3-KMS states w.r.t. the free dynamics (and 
f3 fixed), each characterized by a probability measure n as in (21). Fix such 
a state a;^. It is easy to see that each of the extremal equilibrium states 
a;| has the property of return to equilibrium, i.e., limt^oo '-^p{B*at{A)B) = 

uI{B*B)lj^{A), for all ^, A, B. It follows that 

lim uji'{B*at{A)B) = [ dn{C) ivi{B*B)ui{A). (23) 

The r.h.s. of (23) is in general not equal to uj'^lA): the limiting state depends 
on the initial condition (i.e., on u!p{B*B)). Relation (23) motivates the fol- 
lowing abstract 

Definitions. 1. Let a; be a state on a C*algebra 21, invariant w.r.t. a 
*automorphism group at of 2t. We say that u is asymptotically stable (w.r.t. 

at) if limt^oo uj{B*at{A)B) = cu{B*B)cu{A), for any A, B E ^. 

2. Let ui^, C, E X (a measurable space), be a measurable family of states 
on a C*algebra 21 (in the sense that ^ h- > u!^{A) is measurable for all A G 21) 
and let at be a *automorphism group of 21. Given any probability measure 
II on X we define the state 

f d/iiO 00^. (24) 
Jx 
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We say that the family lu^ is asymptotically stable (w.r.t. at) if, for any 
II, A, B, we have 

lim uj''{B*at{A)B) = / dfx{^) uj^{B*B) uj^{A). (25) 

t^oo 

3. If in 1. is a at)-KMS state then we say uj has the poperty of 
Return to Equilibrium. Similarly, if the uj^ in 2. are at)-KMS states (then 
so is Lu'^) we say the family cu^ has the property of Return to Equilibrium. 

Remarks. 1. More generally one could consider in (24) the case where ji 
is a measure on the space of all states on 21. The present setup is sufficient 
for our purposes. 

2. If oj^ is asymptotically stable w.r.t. at, for all ^, then the family a;^ is 
asymptotically stable. 

3. If B satisfies u^{B*B) = 1 for all ^ (say if B is unitary) then (25) 
means that uj^ is asymptotically stable. In general, the effect of the initial 
condition on the limit state (25) is a redistribution of the relative weights. 

In the above definitions the dynamics of the system is given by a (not nec- 
essarily norm continuous) *automorphism group at of a C* algebra 21. While 
this description applies to free Fermionic or Bosonic heat reservoirs it does 
not in our case of interest, where a Bosonic reservoir is coupled to a small 
system. The problem is that one docs not know how to define the dynamics 
for the coupled system as a *automorphism group of the C*algebra of observ- 
ables (unless the algebra is changed, see [FMl]). One circumvents this issue 
by defining the interacting dynamics, via a converging perturbation series, 
as a *automorphism group of the von Neumann algebra associated with a 
reference state (the uncoupled equilibrium state), sec Section 2.1.2. We shall 
therefore adapt the above definitions to a setting where the dynamics is not 
defined on the level of the C* algebra of observables, but is rather expressed 
as a ( "Schrodinger" ) dynamics of states. 

Definitions. 1'. Let c<j be a state on a C*algebra 21 and denote by 
{'H^,^l^J,^^^) its GNS representation, uj{A) = {fl^^,ni^{A)fli^) . Suppose ctj is 
a *automorphism group of the von Neumann algebra 71^^(21)". We say that 
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LU is asymptotically stable (w.r.t. at) if 

lim {^l^,7r^iB*)ati7r^{A))7r^iB)n^) = u;iB*B)u{A), (26) 

t— >oo 

for all A,B eOl. 

2'. Let a;^, ^ E X {a. measurable space), be a measurable family of states 
on a C*algebra 21 and denote their GNS representations by tt^, f2^). 
Suppose that, for each ^, is a * automorphism group of the von Neumann 
algebra 7r^(2l)". We say that the family cu^ is asymptotically stable (w.r.t. al ) 
if, for any Ijl,A,B e 2t, we have 

^hmy" di^iO {n^,TT^{B*)al{7r^{A))7r^{B)n^) = J di,{i) uj^{B* B) uj^{A), 

(27) 

where fi is an arbitrary probability measure on X. 

3'. If uj in 1'. is a (/9, crt)-KMS state of t^uj{^)" then we say lo has the 
poperty of Return to Equilibrium. Similarly, if the in 2'. are (/3, (t|)-KMS 
states of 7r^(2t)" we say the family uj^ has the property of Return to Equilib- 
rium. 

Remark. In case al{Ti£^{A)) = 7r^{at{A)) for some *automorphism group 
at of 21, the second set of definitions reduces to the first one. 

In this paper we show a (weak) version of relation (27), for equilibrium 
states of the Bose gas with a condensate interacting with a small system, 
where the time limit is taken in the ergodic mean sense and is followed by 
the hmit of small coupling constant, see (40), (39). 

1.2 A condensate coupled to a quantum dot (quantum 
tweezers) and its weak coupling stability 

In this section we outline the model and the stability result. For technical 
detail we refer the reader to Section 2. The small system with which the 
supercritical Bose gas interacts can trap finitely many Bosons - we call it 
therefore a quantum dot. One can imagine the use of such a trap to re- 
move single (uncharged) particles from a reservoir, hence the name quantum 
tweezers (see e.g. [DWRN]). 
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The pure states of the small system are given by normalized vectors 
in C^. We interpret [1,0,..., 0] as the ground state (or "vacuum state"), 
[0, 1, 0, . . . , 0] as the first excited state, e.t.c. The Hamiltonian is given by 
the diagonal matrix 



Hi = diag(0, l,2,...,rf- 1). 



(28) 



Our method applies to any selfadjoint diagonal matrix with non-degenerate 
spectrum. We introduce the raising and lowering operators, and G_, 



G 



+ 





1 

••• 












(29) 



(G+ has ones on its subdiagonal) which satisfy HiG± = G±{Hi ±1). The 
action of G^ (C-) increases (decreases) the excitation level by one. The 
dynamics of an observable A £ B{C^) (bounded operators on C^) is given by 

A ^ a\{A) = e'*^^Ae-'^^\ 

The observable algebra of the combined system is the C*-algebra 



2l = B(C'^)®2II(P), 



(30) 



and the non-interacting dynamics is the *automorphism group of 21 given by 

ai, (31) 



where we now denote the free field dynamics (2) by a^. 

Let /i be a fixed probability measure on [pcrit; oo) x ^-nd consider the 
(/?, q;^)-KMS state 

^Ko- I c^MO^i,/?®^?, (32) 

where oji^p is the (/3, q;^)-KMS state (Gibbs state) of the small system, and a;^ 
is a {I3,a\)-KMS> state with fixed density and phase (see Subsection 1.1.1). 
The subindex "0" in (32) indicates the absence of an interaction. Let H 
denote the (GNS-) Hilbert space of state vectors obtained from the algebra 
21, (30), and the state (32). Furthermore, let VL^^^ G Ti denote the cyclic 
vector in Ti representing the state a;^°o ' tt be the GNS representation 
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of 21 on Ti. Since c<j^°q is invariant under (see (31)) there is a selfadjoint 
operator Cq acting on 7Y, called the thermal Hamiltonian or Liouvillian, 
satisfying 

7r(a*(A)) = e'*^»7r(A)e-**^«, (33) 
£o^^^:S = 0, (34) 

for all A G 21. In order to describe interactions between the small system 
and the Bose gas one replaces the (non-interacting) Liouvillian Cq by the 
(interacting) Liouvillian Cx, which is the selfadjoint operator on H given by 

Cx = Co + AX, (35) 

where A G M is a coupling constant and I is the operator on Ti determined 
by the formal interaction term 

X{G+^a{g) + G_^a*{g)), (36) 

where G± are the raising and lowering operators, (29), and a'^^g) are cre- 
ation (t^ = *) and annihilation operators of the heat bath, smeared out with 
a function g E V, called a form factor. The operator ® a{g) destroys a 
Boson and traps it in the quantum dot (whose excitation level is increased 
by one) and similarly, the effect of G_ (g) a*{g) is to release a Boson from 
the quantum dot. The total number of particles, measured by the observ- 
able Hi + J^^ a*{k)a{k)d^k, is preserved by the interaction (meaning that 
(36) commutes with this operator). Since the quantum dot can absorb only 
finitely many Bosons, the interacting equilibrium state is a (local) perturba- 
tion of the non-interacting one. A physically different situation occurs when 
the condensate is coupled to another reservoir. We would then expect that 
time-asymptotic states are of non- equilibrium stationary nature, a situation 
we investigate in a different work. 

Of course, (36) has a meaning only in a regular representation of the 
Weyl algebra, e.g. the representation tt above, see Subsection 2.1.2. The 
interaction I has the property that the dynamics generated by C\ defines 
a ^automorphism group <t\ of the von Neumann algebra 971^°° C <B(7i) ob- 
tained by taking the weak closure of the algebra 7r(2l). One can show that 
there exists a vector Q'^^l G H defining a (/?, a\)-KMS state on mt^°". We 
call the perturbed KMS state, it satisfies 

^x^Z = 0. (37) 
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As can be seen from (32) the Hilbert space Ti has a direct integral decompo- 
sition, 



I dfiiOH^, (38) 



where Ti^ is the Hilbert space associated with the state cui^p <S> uj^. Similarly, 
the representation vr, the cyclic vectors ^'j^Q, ^Jf") ^^e von Neumann algebra 
971^°'^ and the Liouville operators Cq, C\ are direct integrals over correspond- 
ing fibers (labelled by ^ and weighted with the given measure /x). 

We are thus in the setting of Definition 2' with the following identifica- 
tions: Q.^ is the fiber of tt^ is the fiber of tt and is the *automorphism 
group e**-'^^'«(-)e~**^^'« of 7rg(2l)", where L^,^ is the fiber of C\. Our weak cou- 
pling result on Return of Equilibrium reads: for all /x. A, B 

hm hm 1 fdtu:^{B*ai{A)B)= I dtx{i)u:l,{B*B)u:%,{A), (39) 

where u;| g is the state determined by fi| q, and a\ is the interacting dynamics 
generated by Cx. The expression <j\{A) has to be understood cum grano salis, 
in the sense of Definition 2' (see also (108)). 

The result (39) follows from our Theorem 2.1: for e > 0, ^ G [Pcrit, oo) xS"-*^, 
S e 21 fixed, there is a Aq > s.t. if < |A| < Aq then 



^lim ^J^^4,,{B*al^{A)B)dt - 4^,{B* B)ujI,{A) 



<e\\A\\, (40) 



for all A G 21, and where ^^ is the vector state given by x- expect 
that Return to Equilibrium should hold in the sense of Definition 3' but can- 
not prove it due to reasons we outline in Subsection 2.2.1. 

The proof of Theorem 2.1 is based on our Theorem 2.3 which describes 
the structure kerL^,^, which in turn follows from a positive commutator ar- 
gument combined with a new Virial Theorem, Theorem 2.2. 

We prove (40) under a condition of regularity and "effectiveness" of the 
interaction. Let us close this section by discussing the physical meaning of 
the latter condition. Consider first the Bose gas at critical density Pcrit(/3) 
for some fixed temperature 1/(3 (so that there is no condensate, po = 0). 
Heuristically, the probability of trapping a Boson in a state / in the quantum 
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dot is given by 

' 'G+®a{f)(p®n,e-''"^ip(»n) , (41) 



where (/? an eigenstate of the quantum dot Hamiltonian and the Bose gas is in 
the equihbrium state Q (for the calculation, we put the system in a box and 
r2 is a vector in Fock space with Bosons distributed according to a discrete 
distribution approaching the Planck distribution as the box size increases). 
The interacting Hamiltonian is Hx — Hq-\- X{G+ ® a{g) -\-G-® The 
second order contribution in A to (41), for large values of t, is 

p^ = c ^^,J_^y \mm\ (42) 

where we assume that f{r),g{r) are radially symmetric, and where C is a 
constant independent of f3, f, g. P2 gives the probability of the second order 
process where a Boson gets trapped in the quantum dot; the excitation energy 
is 1 (the quantum dot Hamiltonian (28) has equidistant eigenvalues) and the 
probability density of finding a Boson with energy uj{l) — 1 per unit volume 
is oc (e^ — 1)~^, according to (6). In order not to suppress this trapping 
process at second order in the coupling constant we assume that g{l) ^ 
( "effective coupling" ) . 

Next let us investigate the influence of the condensate. For this we fix 
a density po of the Bose gas and consider very low temperatures {(3 — > 00), 
so that most particles are in the condensate. If Q denotes the corresponding 
state of the Bose gas then we calculate the second order in A of (41) to be 

Q2{t) = C(l - cost)AVol/(0)y(0)|l (43) 

We see from (43) that if g{Qi) — then there is no coupling to the modes of 

the condensate: a physically trivial situation where the condensate evolves 
freely and the small system coupled to the "excited modes" undergoes return 
to equilibrium. In this paper we develop a theory which includes the case 
5'(0) 7^ 0, a situation which could not be handled by approaches developped 
so far. 

2 Definition of the model and main results 

In Section 2.1 we introduce the class of systems considered in this paper 
and we explain the central decomposition of the equilibrium state with a 
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condensate (references we find useful for this are [AW], [LP], [C] and also 
[H]). Our main results are presented in Section 2.2, at the end of which we 
also give the quite short proof of the stabihty result, Theorem 2.1. 

2.1 Definition of model 

We introduce the uncoupled system in Subsection 2.1.1 and present its Hilbert 
space (GNS) description (see (71), (72)) including the uncoupled standard 
Liouvillian Lq, (80). The interaction is defined by an interacting standard 
Liouvillian introduced in Subsection 2.1.2, (96). 

2.1.1 Non-interacting system 

The states of the small system arc determined by density matrices p acting 
on C^. A density matrix is a positive trace-class operator, normalized as 
trp = 1, and the corresponding state 

Up{A) = tr {pA), A e E{C^) (44) 

is a normahzed positive hnear functional on the C*-algebra B{C^) of all 
bounded operators on C^, which we call the algebra of observables. The 
(Heisenberg-) dynamics of the small system is given by the group of *auto- 
morphisms of B{C^) generated by the Hamiltonian Hi given in (28), 

a\{A)^^'"^Ae-''^\ t e R. (45) 

Denote the normalized eigenvector of Hi corresponding to the eigenvalue 
Ej by ipj. Given any inverse temperature < /5 < oo the Gihhs state uji^p 
is the unique /3-KMS state on B{C^) associated to the dynamics (45). It 
corresponds to the density matrix 

Let p be any density matrix of rank d (equivalently, p > 0) and let {ip j}'^~l 
be an orthonormal basis of eigenvectors of p, corresponding to eigenvalues 
< pj < 1, Y^jPj = 1- The GNS representation of the pair (;B(C'^), Wp) 



16 



is given by (Tii, tti, where the Hilbert space 7ii and the cydic (and 
separating) vector fli are 

Hi = C^^C^ (47) 
Qi = ipj0 ifj e (g) C"', (48) 

3 

and the representation map tti : B{C^) B{Hi) is 

7ri(^) =A®1. (49) 
We introduce the von Neumann algebra 

=8(^^)0 led C BiUi). (50) 
The modular conjugation operator Ji associated to the pair (SDTi, fti) is given 

by 

Jl"?/'^ ® 'i/'r = Ci-?/'^ (g) Ci'^^, (51) 

where Ci is the antilinear involution Ci ^jZjipj = '^j^^pj (complex conju- 
gate). According to (48) and (46) the vector fli^p representing the Gibbs 
state LUi^is is given by 

n,,^ = J2 e-^'^^/Vi ® 'Pj e ^1- (52) 

Vtre-^-^i J 

We now turn to the description of heat bath. Its algebra of observables 
is the Weyl algebra 2IJ(X>) over some linear subspace of test functions T> C 

L'^ (R^ , k) . The elements of I) represent the wave functions of a single 
quantum particle of the heat bath. The choice of V depends on the physics 
one wants to describe - in particular, it is not the same for a system of 
Bosons with and without a condensate, as we will see shortly. For fixed T>, 
W{V) is the C*-algebra generated by elements W{f), f eV, called the Weyl 
operators, which satisfy the CCR (1). The ^operation of is given by 

W{f)* = W{—f). The dynamics of the heat bath is described by the group 
of *automorphisms of W{V) 

al{W{f))^W{e''^f), (53) 
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where /i is a selfadjoint operator on L^(M.^,d^k). In the present paper, we 
choose h to be the operator of multiphcation by the function u{k), see (3). 
Our methods can be modified to accomodate for other dispersion relations 
than (3). 

Wc choose the test function space V to consist of all functions / G 
L^{M:^,dH:) s.t. ^;^^(/) exists (see (11)). If po = the r.h.s. of (11) re- 
duces to the product of the first two exponentials (one may then extend 
V to L2(R3^ (1 + p)d^k)), and if in addition /3 ^ oo p{k) = 0) then 
E{f) — e~4 11/11^ is just the Fock generating functional corresponding to the 
zero temperature equilibrium state (in this case one may extend the test 
function space to all of L^(]R^, rf'^/c)). 

The GNS representation of the pair (W{V), 002,13), where uj2,i3 is the equi- 
librium state with supercritical density p (determined by E'^^^{f), (11)), has 
been given in [LP] as the triple (7^2, '^2, ^2), where the representation Hilbert 
space is 

n2 = J^0j^0L\R'^,dpf3^p), (54) 

where = J^[L^(M:\d,^k)) is the Bosonic Fock space over L^(M^,d^/c), and 
the probability measure dpp^-p is given in (19). The cyclic vector is 

Q2 = QjF^njF® 1 (55) 

where is the vacuum in !F and 1 is the normalized constant function in 
L^(M^, dpp^p). The representation map 7:2 : W{T>) — > B{'H2) is given by 

MWif)) = WAVlTpf) (8) W:r{Vpf) ® e'^^^^'^, (56) 

where Wyr = e'-f^'^^) is a Weyl operator in Fock representation, the field 
operator </?^(/) is 

M/) = ^(«^*(/) + M/)) (57) 

and ajF*(/), ajc-(/) are the smeared out creation, annihilation operators sat- 
isfying the commutation relations 

\aT{!\aj:\g)\ = (/,^) , [a^(/), 0^(5)] = [a^*(/), a^*(^)] = 0. (58) 

Our convention is that / 1— > a:F(/) is an antihnear map. The phase $ e R is 
given by 

$(/,{) = (27r)-'VV2(r - p^rit) ((Rc/(0)) cos^^ + (Im/(0)) sm^^) , (59) 
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for ^ = {r,6) G [pcritjOo) x S^. In the absence of a condensate the third 
factor in (54), (55) and (56) disappears and the representation reduces to 
the "Araki- Woods representation" in the form it has appeared in a variety 
of recent papers. We denote this representation by ttq. More precisely, let 
TTo : W{'D) E{T ® T) be the representation 

vro(l^(/)) = Wr{^/lVpf) ® Wr{^J). (60) 

It is well known that the vector 

Qq^Qj.® fijr (61) 

is cyclic and separating for the von Neumann algebra 

9Jto := 7ro(2U(P))" ^ B{T ® T) (62) 

(it defines a KMS state w.r.t. the free field dynamics). We also introduce 
the von Neumann algebra 

m2 := 7r2(2IJ(P))" = mio ® M C B(7^2), (63) 

where }A is the abelian von Neumann algebra of all multiplication operators 
on L^(R^, d//^_p). It satisfies M' — M. The equality in (63) follows from 
this: clearly we have SDTq' ® A4 C VJt^' . so taking the commutant gives DXIq <® 
A4 D 9JI2. The reverse inclusion is obtained from Ij^^j^ <S) C. and 
9Jlo ® ]li2(R2) c 9JI2 (sec [LP]). 

It is well known that the von Neumann algebra DJIq, (62), is a factor. 
That means its center is trivial, 3(SDTo) = ^0 H SDTo' = C. However, we have 
3{Wl2) = {Mo ®M)r) {DJlo ® M), i.e. 

3{m) = tr^r®M, (64) 

so the von Neumann algebra 97l2 is not a factor. This is clear on physical 
grounds since the system has long-range correlations, see Section 1.1.1. One 
can decompose 9Jt2 into a direct integral of factors. The Hilbert space (54) 
is the direct integral 

^2= / c?/^/3,p(0 (65) 
Jr2 
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and the formula (see (55), (56), (60), (61)) 



(66) 

shows that 7r2 is decomposed as 

/■e 

7r2 = / rf/i/3,p(0 TT^, (67) 

where tt^ : W{T>) — > B(J^ ^) is the representation defined by 

TTdWifj) = e-*(^'«)7ro(iy(/)). (68) 

For each fixed ^, 

7r^(2H(P))" = 9Jto (69) 
is a factor. Accordingly we have 



Mo 



r dii0-,{i)mo. (70) 



The GNS representation of (21, cu^^q) (see (32)) is just (7i,7r, fi), where 

TT = TTi ® 7r2 (71) 

n^"^ = ni,;3®02. (72) 

The free dynamics is given by the group of * automorphisms a^, (31). Let 

9Jtr — 7r(2l)" = 9Jli ® 9Jl2 = / c?Af/3,p(0 ^i^^o C B(7i) (73) 

be the von Neumann algebra obtained by taking the weak closure of all 
observables of the combined system, when represented on Ti. To see how 
we can implement the uncoupled dynamics in 7i we use that (for all t G M) 
$(e''^*/,^) = $(/,0, which follows from a;(0) = 0, see (59) and (3). Thus 

^2{al(W(m = r d^i^-^iOe-'^^f'^^nome'^'f)). (74) 
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It is well known and easy to verify that for ^4 e 21, 

(tti (8) 7ro)(a*(^)) = e^*^°(7ri 7ro)(>l)e-^*^°, (75) 

where the self adjoint operator Lq on T^i ^ is given by 

Lo = L1 + L2, (76) 
Li = ifi led - led <S> Hi, (77) 
L2 = dr(u;)®]l^-]l^(8)dr(a;). (78) 

Here dr(a;) is the second quantization of the operator of multiplication by ou 
on L'^(M.^,cPk). We will omit trivial factors 1 or indices c^, t whenever we 
have the reasonable hope that no confusion can arise (e.g. Li really means 
Li (g) lljr (g) 11 jr). It follows from (73)-(78) that the uncoupled dynamics al is 
unitarily implemented in Ti. by 

n{al{A)) = e^*^°7r(A)e-^*^°, (79) 

where the standard, non-interacting Liouvillian Cq is the selfadjoint operator 
on H. with constant (^-independent) fiber Lq, 

>Co = / diip-,{i) Lo. (80) 

The r.h.s. of (79) extends to a *automorphism group ctq of which is 
reduced by the decomposition (73). We write 

4 = / dfip-,{0 al^, (81) 

where cTq^ is the ^automorphism group of 9Jti 9Jlo generated by Lq. As is 
well known, 

^p,o = ® ^0 (82) 

is a (/5, (To^^)-KMS state of 9Jli OJIq. The modular conjugation operator J 
associated to (QJlo, <8) f^o) is 

J = Ji ® Jo, (83) 
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where Ji is given by (51) and where the action of Jq on ^(g)^ is determined 
by antihnearly extending the relation 



JoMWifmo = W:p{^f) (8) W^{y/lT^f)no. (84) 

Jo defines an antihnear representation of the Weyl algebra according to 
W{f) H- > Jo7ro(VF(/)) Jo, which commutes with the representation ttq given 
in (60). We view this as a consequence of the Tomita-Takesaki theory which 
asserts that TIq = JoOJIqJo. 

It follows from (72), (73), (81) that 

^^^S = / df,0-,iO ^0,0 (85) 

is a (/3, (To)-KMS state on 9Jl^°°, and that the modular conjugation operator 
J associated to (97t^°°, Q^°o) is given by 

J = / d///3,p(0 Ji Jo- (86) 

The standard Liouvillian Cq, (80), satisfies the relation 

JCo = -Co J. (87) 

One can choose different generators to implement the dynamics on Ti 
(by adding to the standard Co any selfadjoint element affiliated with the 
commutant (071^°°)'). The choice (80) is compatible with the symmetry 
9Jl^°° = {dJVp"^y, in that it also implements al ^or the antilinear representa- 
tion Jt^{j)J . Another way to say this is that the standard Liouvillian (80) 
is the only generator which implements the non-interacting dynamics 
satisfies 

Co^^^^l = 0, (88) 

see e.g. [BR, DJP]. 

2.1.2 Interacting system 

We define the coupled dynamics, i.e. the interaction between the small sys- 
tem and the Bose gas, by specifying a *automorphism group a\ of the von 
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Neumann algebra (the "perturbed" or "interacting dynamics"). One 

may argue that a conceptually more satisfying way is to introduce a repre- 
sentation independent regularized dynamics as a *automorphism group of 21 
and then to remove the regularization once the dynamics is represented on 
a Hilbert space. This procedure can be implemented by following the argu- 
ments of [FMl], where it has been carried out for the Bose gas without con- 
densate. The resulting dynamics is of course the same for both approaches. 
For a technically more detailed exposition of the following construction we 
refer the reader to [FMl]. 

The interaction between the two subsystems is given formally by (36), 
which we understand as an operator in a regular representation of the Weyl 
algebra, so that the creation and annihilation operators are well defined. We 
could treat interactions which are sums over finitely many terms of the form 
(36), simply at the expense of more complicated notation. 

The field operator (/?(/) — \dt\t=o'^(VV{tf)) in the representation tt, (71), 
is easily calculated to be 

^(/) = / dfi^rpiOMf)^ (89) 

cp^if) = MyrT7/)®]i+]i®MVp7)-*(/>e), (90) 

where $(/, is given in (59), and where <fj:-{f) is given in (57). Define the 
interaction operator by 



V ^ G+0lcd0 (^a^(Vl + pg) 1^+1^0 ar{\fp9) 

-(27r)-3/V2(r-pcrit) W) e*') + adjoint, (91) 

which corresponds formally to 7r(G+ ® a(fi')) + adjoint (apply (71) to (36)). 
V is an unbounded selfadjoint operator on Ti. which is affihated with 
For i e K, A e 9Jl^°° we set 



al{A) = ^{i>^T dh... dtn[e''"^°Ve-''-^°,[--- 

n>0 

[e^^'^°Ve-'^'^°, A] ■■■]]. (92) 



The series converges in the strong sense on a dense set of vectors, for any 
A e !m^°°, A,t e M (see e.g. [FMl]). Since V is affihated with fm^™^ and 
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git/Zo .g «t^o leaves OJt^°" invariant, the integrand in (92) does not change if we 
add to each e'^^'^'>Ve-'^^^° a term -Je^h^^oy^-itjCoj ^ _^it,Cojyj^~itjCo 
(because this term is affihated with the commutant (QJl^""^)'). In other words, 
V in (92) can be replaced hy V - JVJ. The r.h.s. of (92) is then identified 
as the Dyson series expansion of 

gitA^g-itA^ (93) 
where the standard, interacting Liouvilhan Cx is the selfadjoint operator 

Cx^Co + X{V - JVJ) = £0 + AX. (94) 

Subtracting the term JVJ serves to preserve the symmetry (87) when 
adding the perturbation, i.e., we have JCx — —CxJ- It is not hard to 
verify that (93) defines a *automorphism group 

al{A) = e^*^Me-^*^^ (95) 

of 9Jl^°'^ which defines the interacting dynamics. The Liouvilhan Cx is reduced 
by the direct integral decomposition, 

Cx = r dixp-p{i) Lx,^, (96) 
Lx,^ = Lo + XI^, (97) 
where Lq is given in (76) and where we define 

h = I + K^^ (98) 
I = G+ (g) % (8) |a^(Vl + pg) ® 1^ + 1^ ® a:r*{^/pg)j + adj. 

® CiG+Ci ® |a^*(VP^) ® 11^ + 11^ ® a^(Vl + P^)} + adj.(99) 
= (g) led (g) 1^^^ - led «) Ciis:|Ci (g) 1^®^ (100) 
X| = -2(27r)-=^/2VF^^(G'+^e^^ + G_^(0)e-^) (101) 

with Ci given in (51), and where the creation and annihilation operators a^F*, 
ajr are defined by (58). It is convenient to write (compare with (81)) 
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where a\ ^ is the *automorphism group on ® 9Jlo generated by La,^, (97). 

To the interacting dynamics (95) corresponds a /3-KMS state on 
the equihbrium state of the interacting system. It is given by the vector 

^7^- r di^P,p(0 ^$,x, (103) 

where 

= (4,A)"'e-''^^°+^'^'^)/'^^/3,o &n^®T®T, (104) 

with a normahzation factor ensuring that llf^l^ll ~ 1' where is 
obtained by dropping the second term (the one coming with a minus sign) 
both in (99) and in (100). The fact that ^p^i (82), is in the domain of the 
unbounded operator e~'''(^°+'^^f'^)/^, provided 

||^/V'^||l2(m3) < oo, (105) 

can be seen by expanding the exponential in a Dyson series and verifying 
that the series apphed to Q^^o converges, see e.g. [BFS]. It then follows from 

the generalization of Araki's perturbation theory of KMS states, given in 
[DJP], that is a (/3, (j^^^)-KMS state on Mi ® Mq, and that 

La,^(]|;, = 0. (106) 
We conclude that Vtf^ is a (/?, (7^)-KMS state on aJl^^"^, and that Cx^t"^^ = 0. 

2.2 Main results 

We make two assumptions on the form factor g determining the interaction 
(see (36), (91)). 

(Al) Regularity. The form factor is a function in C^(M^) and satisfies 

11(1 + l/^/uJ){k ■ Vfe)^5f||L2(]R3,d3fe) < oo, 

for j = 0, . . . , 4, and || (1 + ufgWL^w^d^k) < oo. 

(A2) Effective coupling. We assume that Jg2da \g{l,a)\'^ ^ 0. Here, g is 
represented in polar coordinates. 
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Remarks. 1) Condition (Al) is used in the application of the virial theorem 
- we choose the generator of dilations ^{k ■'Vk + Vfc ■ k) to be the conjugate 
operator in the theory. 

2) Condition (A2) is often called the Fermi Golden Rule Condition. It 
guarantees that the processes of absorption and emission of field quanta by 
the small system, which are the origin of the stability of the equilibrium, 
are effective, see the discussion in the introduction. In a situation where the 
spectrum of Hi has gaps between neighbouring eigenvalues, condition 
(A2) is replaced by rninj Jg2da\g{uj~^{Aj),a)\^ ^ 0. This means that the 
form factor should couple the modes of the gas which induce transitions of 
the small system. 

Theorem 2.1 (Weeik coupling return to equilibrium). Assume con- 
ditions (Al) and (A2). Let e > 0, ^ e [pcrit, oo) x 5"^ S e 21 be fixed. There 
is a Ao(e, ^,B) > s.t. ifO< \X\ < Ao(e, ^, B) then 

}^^f£dt (7r,(S)n|,,,4,(7r,(A))7r,(S)0|,,) 

<e\\Al (107) 

for all >1 e 21. The coupled KMS state a;| ^ is determined by the vector (104)- 
Set 

u^{B*ai{A)B):= [ d,,{0 (TidmU^ 4xMA))^dBH,x) , (108) 

where fj, is any probability measure supported on [pcrit^oo) x S^- If follows 
from (107) that 

lim lim l fdtu;>^{B*ai{A)B)^ [ df^iO u;p,o{B*B)u;f,,o{A), (109) 

where cup^o is the uncoupled KMS state determined by the vector ^Ip^, (82). 

Remark. We expect the r.h.s. of (107) to be zero for A sufficiently 
small. This has been proven to hold for systems without a condensate (with 
varying conditions on the interaction and varying modes of convergence) in 
several papers, see [JPl, BPS, Ml, DJ, FM2]. The obstruction to applying 
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the strategies of these papers is that they all need the condition that either 
g{0) = 0, or g{k) ~ as \k\ — >• 0. The first case is uninteresting in 

the presence of a condensate (no couphng to the modes of the condensate!), 
and the second type of form factor does not enter into the description of a 
system with a condensate (see (59)). We refer to Subsection 2.2.1 for a more 
detailed discussion of this point. 

In order to state the virial theorem and to measure regularity of eigen- 
vectors of La,^, (97), we introduce the non-negative selfadjoint operator 

A = dr(c^)(^ V+]l^(8dr(a;), (110) 

where dr(ti;) is the second quantization of the operator of multiplication by 

cl-'(A;) on L'^(M.^,d^k), c.f. (3). The kernel of A is spanned by the vector 
Qq = ® (c.f. (61)) and A has no nonzero eigenvalues. The operator 
A represents the quadratic form i[Lo,/l], the commutator of Lq with the 
conjugate operator 

A = dr(ad)(8 V- V0dr(ad), (111) 
where Od is the selfadjoint generator of dilations on L^(R^, (Pk), 

ad^i(k-Vk + ^y (112) 

The formal relation A = i[Lo,A] follows from i[c(;,ad] = oo (for relativistic 
Bosons, see (3); in the non-relativistic case we have i[uj, ad/2] = w). The self- 
adjoint operator representing the quadratic form i[Lx,^, A] is easily calculated 
to be (see (99)) 

Ci = A + A/i (113) 
Ii = 6*+ ® led ® |ajF(adyiTp5') <8) Ijr - IjF <8) a:r*(adv^^)| + adj. 

® CiG+Ci ® ^a^*{ad^/pg) ® 11^ - 11^ ® a^(adVl + p^)} + adj. 

(114) 

Similar expressions are obtained for the higher commutators of La,^ with 
A, see Section 3. Assumption (Al) guarantees that (1 + l/^/uj){ady ^/pg 
and (1 + l/y/Lj){ady^/TTpg are in L'^{R^,d^k), for j = 0,...,4, so the 
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commutators of La,^ with A, up to order four, are represented by selfadjoint 
operators (satisfying the technical requirements needed in the proof of the 
virial theorem). 

Theorem 2.2 (Virial Theorem, regularity of eigenvectors of Lx^^)- 
Assume condition (Al) and let ^ G [pcrit, oo) x 5*^ he fixed. If ip is an eigen- 
function of L^^^ then ip is in the form domain of Ci, (113), and 

(^,Ci^) = 0. (115) 

There is a constant c which does not depend on ip, ^ nor on {3 > Pq, for any 
/5o > fixed, such that 

||AV2V||<c|A|||V'||. (116) 

Remarks. 1) Relation (115) seems "obvious" from a formal point of 
view, writing Ci = ^[La,^, A] = i[Lx^^ — 6,^4], and using that {Lx^^ — e)* — 
Lx,^ — e, where Lx^^i^ — eip. A proof of (115) is certainly not trivial, though, 
and considerable effort has been spent by many authors to establish "Virial 
Theorems" (see e.g. [ABG] and [GG] for an overview, and also [Ml], [FMl] 
for approaches similar to ours). 

2) The regularity bound (116) follows easily from (115) and (113) and 
from the standard fact that Ii is infinitesimally small relative to A^/^ (Kato), 
so that = (ipjCiip) > (1 — e) {ip,Aip) — — c||-?/'|p, for any e > 0, for some 
constant c independent of C, and j3, as mentioned in the theorem. We refer 
for a more complete exposition of this to [FMl]. 

We prove Theorem 2.2 in Section 3.2 by showing that the hypotheses 
leading to Theorem 3.2, a more general result, are satisfied in the present 
situation. Our next result describes the structure of kerL^,^. Let P(A < x) 
stand for the spectral projection of A onto the interval [0, x]. 

Theorem 2.3 (Structure of the kernel of I^a,^)- Assume Conditions 
(Al), (A2) and let ^ G [pcriti oo) x he fixed. There is a number Aq > s.t. 
if < \X\ < Aq then any normalized ip G ker(LA.^) satisfies 

\\P,,^P{A<\Xm>l-0{X'), (117) 

where P^^p is the projection onto Cfii^^ (see (52)) and 0{}^) is a vector whose 

norm, which is independent ofip, tends to zero in the limit A — (uniformly 
in ^ in any compact set and in jS > I3q, for any /?o > fixed). The constant 
Ao is uniform in ^ in any compact set, and in 13 > (3o, for any fixed (3o > 0. 
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Our proof of this theorem, given in Section 5, reUes on a positive com- 
mutator estimate and Theorem 2.2. Expansion (117) implies that the only 
vector in the kernel of Lx,^ which does not converge weakly to zero, as A — > 0, 
is the interacting KMS state (lO^). This information on the kernel of 
Lx,^ alone enters our proof of Theorem 2.1. 

Corollary 2.4 Assume Conditions (Al) and (A2) and let P^^ the projection 

onto the subspace spanned by the interacting KMS state (104)- Let ^ G 

[pcrit, oo) X 5^ be fixed. Any normalized element ip G ker(LA,j) fl (Ran pI^)"*" 
converges weakly to zero, as A — > 0. The convergence is uniform in ^ in any 
compact set and in (5 > for any Po > fixed. 

We prove the corollary in Section 5. The virial theorem wc present in 
Section 3, Theorem 3.2, is applicable to systems without a condensate, in 
which case one is interested in form factors g which have a singularity at the 
origin. Theorem 3.2 can handle a wide range of such singularities (see the 
remark after Theorem 2.4) and is therefore relevant in the study of return to 
equilibrium and thermal ionization for systems without condensate, as will 
be explained in [FM3]. 

Theorem 2.4 (Improved Virial Theorem for systems without con- 
densate). Let Lx be the Liouvillian of a system without condensate, Lx = 
Lq + \I (i.e., = 0), see (98) and suppose that the form factor g is in 
C^(M^\{0}) and satisfies the condition 

(l + l/V^)(ad)-'VrT7^, {l + l/V^){a,,i)-'^pg E L^R^d^k), (118) 
{l+u;)\aay^/lVpg, {1 + ufia^y g eL\R\d^k), (119) 

for j = 0, ... ,4. Then the conclusions (115), (116) of Theorem 2.2 hold. 

Remark. An admissible infrared behaviour of g satisfying (118), (119) is 
g{k) ~ |A;p, as \k\ ~ 0, with p > —1/2 for relativistic Bosons (c.f. (3)). The 
range of treatable values of p obtained in previous works, [M1,DJ,FM1,FM2], 
is p = —1/2,1/2,3/2, p > 2. Theorem 2.4 fills in the gaps between the 
discrete values of these admissible p. 

The proof of Theorem 2.4 is the same as the one of Theorem 2.2, see 
Section 3.2. 
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Theorem 2.5 Assume the setting of Theorem 2.4, that (A2) holds and that 
\g{k)\ < c\k\P, for \k\ < d , for some constants c,c', and where p > —1/2 (for 
relativistic Bosons, and p > for nonrelativistic ones). There is a number 
Xq > s.t. i/0 < |A| < Ao then any normalized ip e ker{Lx) satisfies 

||Pi,/?P(A<|A|)V||>l-0(A°), (120) 

where 0(A°) is a vector whose norm is independent of ip and vanishes in the 
limit A — >• (uniformly in P > Po, for any (3q > Q fixed). The constant Aq 
does not depend on P > Po, for any fixed Pq > 0. 

We give the proof Theorem 2.5 together with the proof of Theorem 2.3 
in Section 5. 

2.2.1 "Weeik coupling stability" v.s. "return to equilibrium", and 
relation with infrcired regulcirity of the coupling 

A central tool in our analysis of the time- asymptotic behaviour of the system 
is the virial theorem, whose use imposes regularity conditions on the interac- 
tion. In particular, wc must be able to control the commutators of L^,^ with 
the conjugate operator A of degree up to four (sec Section 3.1). Depend- 
ing on the choice of A this will impose more or less restrictive requirements 
on the interaction. A convenient choice for A is obtained by representing 
T ®T T{V-{^ X S'^.du X da)) and choosing A = idr(9„) {translation 
generator). This choice, introduced in [JPl], has proven to be very useful 
and was adopted in [Ml, DJ, FMl, FMS, FM2]. The commutator of the non- 
interacting Liouvillian Lq = dT{u) with A (multiplied by i) is just N = dr(]l), 
the number operator in JF(L^(M x S'^,dux da)), which has a one-dimensional 
kernel and a spectral gap at zero. We may explain the usefulness of the gap 
as follows. If one carries out the proof of Theorem 2.3 with the translation 
generator as the conjugate operator then the role of A, (110), is taken by 
N, and relation (117) is replaced by l|Pi,^P(7V < |A|)^|| > 1 - 0(A°). But 
for |A| < 1, P{N < |A|) = |r2o)(r2o| is just the projection the span of 
(product of two vacua in (8) J^), so one has ■0 = ^/3,o + 0{X^), where Qp^ 
is the non-interacting KMS state. Since Qj^fi is close to n| ^^ for small values 
of A, this means that there are no elements in the kernel of L\^^ which are 
orthogonal to provided |A| is small enough, i.e., the kernel of Lx^e^ has 

dimension one and consequently return to equilibrium holds. 
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The disadvantage of the translation generator is that its use requires 
(too) restricitve infrared regularity on the form factor. Indeed, the j-th com- 
mutator of the interaction with the translation generator involves the j-th 
derivative of the fuction , ^ , so an infrared singular behaviour of this 

function is worsened by each application of the commutator (and we require 
those derivatives to be square integrable!). As a result, the case g{0) ^ 
cannot be treated. 

The remedy is to develop the theory with a conjugate operator A which 
does not affect the infrared behaviour of , ^ in a negative way. The 

choice (111) {dilation generator) is a good candidate (one could as well take 
operators interpolating between the translation and the dilation generator). 
The disadvantage of using the dilation generator is that its commutator with 
the non-interacting Liouvillian gives the operator A, which still has a one- 
dimensional kernel, but does not have a spectral gap at zero. This prevents 
us from showing that the kernel of La,j is simple. We can only prove (117) 
which allows us only to show stability of 0Jp°Jl, in the sense of Theorem 2.1, 
but not return to equilibrium. 

We remark that the dilation generator has been used in [BFSS] to show 
instability of excited eigenvalues in zero-temperature models. We expect 
that it is a relatively easy exercise to modify the techniques of [Ml] and 
show absence of nonzero eigenvalues of Lx,^ (which we view as the "excited 
eigenvalues" in the positive temperature case) by using the dilation instead 
of the translation generator. Notice though that if one succeeds to show that 
the kernel of L^,^ is simple, then one knows automatically that Lx,^ cannot 
have any non-zero eigenvalues, see e.g. [JP2]. 

2.2.2 Proof of Theorem 2.1 

Fix Tj > 0,C, E [pcrit, cxd) X S^,B G 21 and choose an element b'^^^ G 7r^(2l)' 
s.t. Ti^{B)ilpfi — h'^^^Vtisfi = 0{7]) {^p,o is cyclic for the commutant 7r^(2t)'). It 
follows that 7r^(S)Q| - b'^^^Qp^ = 0{r]+\\B\\X°) (we use that - = 
0(A°), see [FM2]) and consequently 

(7r^(B)f^|_„ e^*^^^%^(^)e-*^^^%4(5)0|^,) (121) 
= (7r,(S)Q|,„6^,/*^^.%,(A)Q|,,) \\B\\{r] + i\\B\\ + ll&^JDA'^))- 
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We use that h'^,^ commutes with e**^^'«7r^(74)e~**''^^'«, that L^.^f^l;^ = and 
the above estimates. 

The von Neumann ergodic theorem tells us that the ergodic average, 
^ Jq^, of the first term on the r.h.s. of (121) converges in the limit T — > oo 
to 



oo 

+ J2 Uim$,xAA>) {4xMmix) , (122) 



where IIa,^ is the projection onto the kernel of L^,^ and f2| U {'4'j,x}j>i is an 
orthonormal basis of ker L;^^^. It follows from Corollary 2.4 that the series in 
(122) converges to zero as A — > 0. Moreover we have 

{n^{B)n$^„ = 4AB*B)+o{m {v+m\ + \wa)>^')), (123) 

where we use the estimates given at the beginning of the proof. The com- 
bination of (121), (122), (123) implies that there exists a Ai(^,77) > s.t. if 
< |A| < Ai(^,77) then 



T^oo h [ (^«(^)^1a, e^*^^.%e(^)e-^*^^'^7re(S)Q|^, 



ulAB*B) + R, ujI^{A) + R,, (124) 



where 



R, = 0(\\B\\{r^+m\ + \\b'a\)>^')). 
R2 = 0[\\A\\\\B\\{r, + {\\B\\ + \\bl^\\)X')). 

Given e > (as in the theorem) we can choose first rj small and then A small, 
in such a way that Ri < e/2 and R2 < ||A||e/2. The l.h.s. of (124) minus 
cUp ^{B*B)uj^^ ^{A) is then bounded in absolute value from above by e\\A\\. ■ 
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3 Another abstract Virial Theorem with con- 
crete applications 

In this section we introduce a virial theorem in an abstract setting covering 
the cases of interest in the present paper (but which is general enough to 
allow for future generalizations). The virial theorem developed in [FMl], 
where the dominant part of [L, A] commutes with A, does not apply to the 
present situation; here the leading term of [[L, A], A] is L. 

3.1 The abstract Virial Theorem 

Let 7i be a Hilbert space, V d Ti a, core for a selfadjoint operator F > 1, 
and X a symmetric operator on T>. We say the triple (X, Y, T>) satisfies the 
GJN (Glimm-Jaffe-Nelson) Condition, or that {X,Y,T>) is a GJN-triple, if 
there is a constant A; < oo, s.t. for all t/j e T>: 

WW < k\\Y^\\ (125) 
±i{{Xijj,Yijj) - {Yijj,Xijj)} < k{i^,Yi^). (126) 

Notice that if {Xi,Y,V) and (Xs,!",©) are GJN triples, then so is (Xi + 
X2,Y,V). Since Y >!, inequality (125) is equivalent to 

IIX^II <A;i||r^||+A:2||^||, 

for some /ci, /c2 < oo. Condition (125) is phrased equivalently as "X < /cF, 
in the sense of Kato on I>" . 

Theorem 3.1 (GJN commutator theorem) If {X,Y,'D) satisfies the 

GJN Condition, then X determines a selfadjoint operator ( again denoted by 
X), s.t. V{X) D T^{Y). Moreover, X is essentially selfadjoint on any core 
for Y , and (125) is valid for all ip e 'D{Y). 

Based on the GJN commutator theorem we next describe the setting 
for our general virial theorem. Suppose one is given a selfadjoint operator 
Y > 1 with core V <Z Tt, and operators L,A,A> 0, D, On, n = 0, . . . , 4, all 
symmetric on T>, and being interrelated as 

{ip,Dij) = t{{Lip,A^)-{Aip,L^)} (127) 
Co = L 

{^,Cn^l^) = z{(C„_i<^,^^)-(A(^,C„_i^)}, n=l,...,4, (128) 
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where Lp^ip ^ We assume that 

(VTl) (X,y,r>) satisfies the GJN Condition, for X = L,A,D,Cn. Conse- 
quently, all these operators determine selfadjoint operators (which we 
denote by the same letters). 

(VT2) A is selfadjoint, V C T>{A), e^*^ leaves V{Y) invariant, and 

^tAy^-itA < ^gfc'ltly^ f g (129) 

in the sense of Kato on D, for some constants k, k'. 

(VT3) The operator D satisfies D < kA^/^ in the sense of Kato on T>, for some 
constant k. 

(VT4) Let the operators T4 be defined as follows: for n = 1, 3 set C„ = A + Vn, 
and set C2 — L2 + V2, C4 — L4 + V4. We assume the following relative 
bounds, all understood in the sense of Kato on T>: 

Vn < A;A^/^ forn = 1, . . . , 4, (130) 
L4 < kA, (131) 
L2 < kA\ for some r > 0. (132) 

Remark. The invariance condition e'*'^X>(y) C T>{Y) implies that the bound 
(129) holds in the sense of Kato on T>{Y), see [ABC], Propositions 3.2.2 and 
3.2.5. 

Theorem 3.2 (Virial Theorem) We assume the setting and assumptions 
introduced in this section so far. If ijj & Ti. is an eigenvector of L then ijj is 
in the form domain of Ci and 

(Ci)^ = 0. (133) 
We prove this theorem in Section 4. 
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3.2 The concrete applications 

The proofs of Theorems 2.2 and 2.4 reduce to an identification of the involved 
operators and domains and a subsequent verification of the assumptions of 
Section 3.1. Let us define 

P = ToiC^C^^ d^k)) ® :^o(Co°°(R^ d^k)), (134) 

where is the finite-particle subspace of Fock space. Take 

y = dr(u; + l)(8)V + ]l^(g)dr(a; + l) + ]l, (135) 

and let the operators L,A,A of Section 3.1 be given, repectively, by the 
operators L^,^ (see (97), or Lx in the case of Theorem 2.4), (110), and (111). 

It is an easy task to calculate the operators Cj; Ci is given in (113), 
6*2 = ^2 + ^-^2! C's = A + AJ3, C4 = L2 + A/4, where L2 is given in (78), and 
where the Ij are obtained similarly to Ji (see (114)). The operator D, (127), 
is just iA[/, A]. It is a routine job to verify that Conditions (VT1)-(VT4) 
hold, with Vji — In and L4 = L2, r = 1. To check Condition (VT2) one can 
use the explicit action of e**^, see also [FMl], Section 8. 

4 Proof of Theorem 3.2 

Before immersing ourselves into the details of the proof we present some facts 
we shall use repeatedly. 

• If a unitary group e^^^ leaves the domain T>{Y) invariant then there 
exist constants k,k' s.t. HFe'^^V^H < ke'''^^\\\Y'ip\\, for all t/' e 1^{Y). 
Moreover, if (X, Y, V) is a GJN triple then the unitary group e'*^ leaves 
ViY) invariant. 

• Let (X, F, V) and (Z, F, V) be GNS triples, and suppose that the 
quadratic form of the commutator of X with Z, multiplied by i, is 
represented by a symmetric operator on D, denoted by i[X,Z], and 
that {i[X, Z],Y, V) is a GNJ triple. Then we have 

^itx^^-itx _ ^ ^ f e'*'^i[X, Z]e-^*i^. (136) 
Jo 
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This equality is understood in the sense of operators on T>{Y). Of 
course, if the higher commutators of X with Z also form GJN triples 
with Y, V then one can iterate formula (136). 

We refer to [FMl] and the references therein for more detail and further 
results of this sort. Let us introduce the cutoff functions 

Mx) = r dye-y\ f{x)^e-y''\ (137) 

J — oo 

9 - gl (138) 

where gi G C^((— 1,1)) satisfies gi{Q) = 1. The derivative (/i)' equals p 
which is strictly positive and the ratio {f'Y/ f decays faster than eponentially 
at infinity. The Gaussian / is the fixed point of the Fourier transform 

/(s) = (27r)-V2 [ dxe-'^^f{x), (139) 

i.e., /(s) = e~*^/^, and we have (/i)' = isfi = p which is a Gaussian itself. 
This means that /i decays like a Gaussian for large |s| and has a singularity 
of type s"^ at the origin. We define cutoff operators, for i/, a > 0, by 

^1,. = g^{vK)^{2T:)-^l^ f dsg,{s)e^^^^ (140) 

7r 

9u = (141) 

/„ = f{aA) = (27r)-^/2 j dsf{s)e'^-^. (142) 

JR 

Since /i has a singularity at the origin, we cut a small interval {—rj^rf) out 
of the real axis, where > 0, and define 

fl^^a-\2n)-'/' f dsUs)e'^'^\ (143) 

where we set = M\{—ri^ri). Standard results about invariance of domains 
show that the cutoff operators g^^ /q, /{'^ are bounded sclfadjoint operators 
leaving the domain T>{Y) invariant, and it is not hard to see that < 
k/a, uniformly in r] (see [FMl]). 
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Suppose that ip is a normalized eigenvector of L with eigenvalue e, Lip = 
eip, \\ip\\ = 1. Let ip E H he s.t. ip = {L + i)^V ^^'^ ^^t {i^n} G V he a. 
sequence approximating (p, ^ (p- Then we have 

V^„=(L + i)"Vn — ^'ip, n^oo, (144) 

and ■0n £ ^(^)- The latter statement holds since the resolvent of L leaves 
'D(Y) invariant (which in turn is true since {L,Y,T>) is a GJN triple). It 
follows that the regularized eigenfunction 

is in V{Y), and that ipa,v,n — > 'i/', as a, z/ — and n ^ oo. It is not hard to 
see that (L — e)'^„ — >• as n — >• cxo, a fact we write as 

(L-e)V'n = o(n). (146) 

Since /^,^ leaves ViY) invariant, and since ViY) C 'D{L), the commutator 
— i[/{|^,L] is defined in the usual (strong) way on T>{Y). We consider its 
expectation value in the state gu'ipn £ ^(^)) 

-'(i/;".ii),„,.=-Hi/;;«.i • ("7) 

The idea is to write (147) on the one hand as {Ci)^^^^ modulo some small 
term for appropriate a, u, n ("positive commutator"), and on the other hand 
to see that (147) itself is small, using the fact that (L — e)'0 = 0. 
The latter is easily seen by first writing 

(L - e)g^ilJn = gu{L - e)'il)n + gi,v[L, gi,u\i^n + [L, gi,u\gi,ui^n (148) 
and then realizing that, due to condition (VT3), 

and similarly, [L, giA9i,u — O {y/i'), so that 
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Next we figure out a lower bound on (147). A repeated application of 
formula (136) gives, in tlie strong sense on V{Y), 



f dsfi{s)e'"^^ [ dsi [ ' ds2 dss [ ' ds4e-'''''^W''''^ 
Jr„ Jo Jo Jo Jo 



(27r)V^ 



+Rr,,iC^ + ^i?^,2C2 + ^Ca, (150) 
where we use that 

(27r)-^/2 / ds {isYf{s)e'''' = f^''\x), 
Jr 

and where we set f[,^ — (/i)'(q;A), = {fi)"{aA), e.t.c, and 

Rr,,n = -z(27r)-V2 T (s)e^^"^. (151) 

J-T) 

Using that q, = ^{aA) — and applying again expansion (136) yields 

n,aCl = faClfa + ^a/a/;.C^2 + ^^/a/^Ca (152) 

~ (2^-^" / /(s)e*^-^ 2' rfsi j^'' ds2 £ dss e--3°^C3e-3°^. 
Plugging this into the r.h.s. of (150) and using that /('q, = '^fafa^ we obtain 
-i{[n.«L]),.,, (153) 



We take the real part on the r.h.s. for free since the l.h.s. is real. The 

error term in (153) is obtained as follows. Clearly we have = O [r]) and 
condition (VT4) gives Cngu = O (z/"'' + z/~^/^) , which accounts for the term 
O {rj/ + 7]/ \fv). The term O [p? j v) is an upper bound for the expectation 
of the terms in (150) and (152) involving the multiple integrals, in the state 
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Qpipn- For instance, the contribution coming from (150) is bounded above as 
follows. Due to condition (VT4) we have 

which gives the following upper bound on the relevant term: 

1~ 



a'^ I ds 

^77 



Ms) 



The integral is finite because fi has Gaussian decay. 

Our next task is to esimtate the real part in (153). It suffices to consider 
the terms 



a'Reia^Cs).^ and a'Re{iQ'Cs) 



(154) 



because f1'^ = 2(/^)^ + 2/^'/a- Let us start with the first term in (154). 
Using the decompostion C3 = A + V3 and the relative bound of V3 given in 
(VT4) we estimate 



We bound the first term on the r.h.s. from above as 



(155) 



a 



Re (/:A/a),,^„ < aW^'OM ||A^/'V'a,.,n|| (156) 



and use that 



to see that for any c > 0, 



O 



a 



a 



Re(/:A/„)^^^^ <- + c{A) 



(157) 
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Choose c = for some ^ > to be determined later. Then, inserting 

again a term Vi into the last expectation value (by adding a correction of 
size O (a^'^^ / , we get 



fl55)| <a'+^ 



{C-. 



1' V'a.i'.n 



+ 0{ — + — + ^ 



3-« 



(158) 



Next we tackle the second term in (154). The Gaussian / is strictly positive, 
so we can write 



a^Re /^A/a\ + O + ^ ) , (159) 



where we have taken into account condition (VT4) in the same way as above. 
It follows that 



1(159)1 < o? 



Ja 



and proceeding as in (156)-(157) we see that 



(Ci) 



^ I a a a' 



Estimates (158) and (160) together with (153) give the bound 



(160) 



-0[^ + — + 



\^/v V 



+ 



+ — + 



(161) 



We combine this upper bound with the lower bound obtained in (149) to 
arrive at 



(l-0(a^+^)) 

^ + o(n) a 



(162) 



= 



a 



— + 



a 



3-? 
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Choose a so small that 1 — (q;^"*"^) > 1/2 and take the limits — > 0, n — > oo 

to get 



O 



+ ^ + — + 



+ 



V 



(163) 



Take for example ^ = 1/2, i/ = i/(q;) = oc'l^. Then the r.h.s. of (163) is 
0(q;V4), so 

hm (Ci) , = 0. 

Since the operator Ci is semibounded its quadratic form is closed, hence it 
follows from fagu{a)i^ — > ■0, q; — > 0, that is in the form domain of Ci and 
that (Ci)^ = 0. ■ 



5 Proofs of Theorem 2.3 and of Corollary 2.4 

In order to alleviate the notation we drop in this section the variable ^ la- 
belling the fiber in the decomposition (65) (imagining ^ e [pcrit, oo) x to 
be fixed). The operator Lx^e^, (97), is thus denoted 

Lx = L^ + X{I + K), (164) 

where / and K are given in (99), (101). In parallel we can imagine that 
K = Q and that Condition (Al) is replaced by (119). 
Let e, p, 6* > be parameters. Set 

Pp = P^P{K<p) (165) 
Po = P{Li^Q) 

Aq = i9\{PpIRl-RllPp) (166) 

Re — PpRe 

R, = (L^ + e2)-V2 (167) 

where Pp = 1 — Pp. We also set Pq = 11 — Pq- The product in (165) is 
understood in the sprit of leaving out trivial factors {Pp — Pq ^ P{A < p)). 
We also define the selfadjoint operator (c.f. (113), (114)) 

P = Ci + i[Lx, Ao] = A + Ii + i[Lx, Ao], (168) 
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where the last commutator is a bounded operator. Let us decompose 



B = PpBPp + PpBPp + 2RePpBPp. (169) 

Our goal is to obtain a lower bound on {B)^^, the expectation value of B 
in the state given by the normalized eigenvector t/j^ of L^. We look at each 
term in (169) separately. In what follows we use the standard form bound 

A/i>-^A-0(A2), (170) 

and the estimates || A^/^-?/'a|| = O (A), \\PoP{A < p)4'\\\ = O (A). The former 
estimate follows from Theorem 2.2 (or Theorem 2.4 for the system without 
condensate) and the latter is easily obtained like this: let x ^ C'o°(^) be such 
that < X < 1, x(0) = 1 and such that x has support in a neighborhood 
of the origin containing no other eigenvalue of Li than zero. Then, for p 
sufficiently smaU, we have PqP^A < p)x{Lq) = 0, so PoP{A < p)^px — 
PoP{A < p){x{Lx) — x{Lq))iJjx = O (A), by standard functional calculus. 
Taking into account (170) we estimate 

(PpBPp)^^ 

> -e\^{Pp[I + K,PpIR]-R]lPp]Pp) -o{\^) 

= 2e\^{PpIR]lPp) +e\^lPpIR]KPp + PpKR]lPp) -O(A') 

> 2eX'{PpIRllPp)^^-^-^0[^ + e), (171) 

where we use in the last step that Pp — PqP{A < p) +P(A > p) to arrive at 
\\PpIR]KPp\\ = \\PpIR^;PoP{A < p)KPp\\ < c. 

The last estimate is due to \\RfPoP{A < p)|| < c and \\PpIP{A < p)\\ < c 
(this follows in a standard way assuming condition (105)). 
Next we estimate 

(^P^^<^)^x ^ \ (^P^)^x - 2^-^'^^ + ^WRI)^^ - O (A^) (172) 
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and 

(Pp{I + K)P,IrI)^^ = \\PM\' O (^^\\P,IR,\\^ = O (^-^^ , (173) 

where we use \\Pp^x\\ < \\PoP{A < p)M + II ^(A > p)V'a|| = O (A/^), 
and ||Pp/i?e|| = O {l/^/e). The former estimate follows from the observations 
after (170) and from ||P(A > p)^a|| < l/y/p\\P(A > p)A^/^t/^All = 0(X/^). 
The estimate ||Pp/i?£|| = 0(1/ y/e) is standard in this business, it follows 
from PpIR^IPp — O (1/e) (see e.g. [BFSS] and also the explanations before 
(179) here below). Combining (172) and (173), and taking into account that 
{PA)^, > im > P)A>^, > p (P(A > p)}^^ > p({P,)^^ - O (A^)) gives 

iP,BP,)^^> I (P,)^y-^0 + (174) 

Our next task it to estimate 

{P,BPp)^^ = A (Pp/iP,)^^ - ^A (Pp{L,PpIRl - IR%)Pp)^^ . (175) 

It is not difficult to see that 

<W.>^, = (P,/iPoP(A<p)>^^ + (P,/iP(A>p))^^ 

= O(A) + O(||(70aA-^/ll|A^/^V'A||) 
= 0(A), 

where (/i)a nieans that we take in Ii only the terms containing annihilation 
operators (see (99)) and where we use ||(/i)aA~^/^|| < c. The second term on 
the r.h.s. of (175) is somewhat more difficult to estimate. We have 

^A (Pp{LxPpIRl - IR%)Pp) 

= -0X' ((/ + K)PpIR])^^ - ex (PpIRlLoPp)^^ 

+eX^(Pp{{I + K)PpIRl-IRlPp{I + K))Pp) , (176) 

where the first term on the r.h.s. comes from the contribution IPpLqII\) 
in the l.h.s. by using that PpL^ = LgPp = LxPp — X{I + K)Pp and that 
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Lxil^x = 0. We treat the first term on the r.h.s. of (176) as 



((i + k)p,irI] 

= ((/ + K)P,IR^PoP{A < p))^^ + ((/ + K)P,IRlP{X > p))^^ 

= O (A||PoP(A < + O (e-^||(/i)aA-'/1 ||A'/Va||) 

= o(a + A^. (177) 

The second term on the r.h.s. of (176) has the bound 

{PpIRlLoP,)^^ = {P,IRlL,P,P{X<p))^^ + {P,IRlL,P{K>p))^^ 
= O (A) + O (||P,(/),i2,A-V2p(A > p)\\ ||AVV,||) 

= O (^) , (178) 
where we use that (with (7)^ = ((7)a)*) 

||Pp(/)aP.A-V2p(A > p)||2 ^ ||p(A > p)A-V2p,(/),p^||2 = O Q . 

The latter bound can be shown by using the explicit form of the interaction 
/, given in (99), and by using standard pull-through formulae to see that a 
typical contraction term in Pp(7)ai?gA~^P(A > p)(I)cPp has the form 

and is thus bounded from above, in norm, by a constant times 1/e, provided 
p > —1/2 (recall that p characterizes the infrared behaviour of the form 
factor, see Theorem 2.5; in the case of the system with condensate we have 
p — 0). To see this use (A + \k\)~^ < \k\~^, and then standard estimates 
which show that the resulting operator is of order €~^; the mechanism is 
that the main part comes from the restriction of the operator to Ran PoPn,, 
(p — 0) and there the resolvent, when multiplied by e, converges to the Dirac 
delta distribution 5{Li ±a;), so the integral is 1/e times a bounded operator. 
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See also Lemma 6.4 of [BFSS]. 

Next we estimate the third term in the r.h.s. of (176) as 



Pp((7 + K)PpIRl - IR^I + K))P 

= O (e-3/2||P^^,||) + O {\\PpIRl{I + K)Ppij, 

= ° (^) • 

where we use again that ||Pp/it!e|| = 0{l/y/e), WPp-^xW — 0{\/^/p), and 
that ||PpJi?^J|| = O (1/e^). Collecting the effort we put into estimates (177), 
(178) and (179) rewards us with the bound 

(W.>-^0(l.^+jA^), (180) 

which we combine with (171) and (174) to obtain 
Wfc (181) 



■T^2 . 



The non-negative operator PpIR^IPp has appeared in various guises in many 
previous papers on the subject ("level shift operator"). The following result 
follows from a rather straightforward calculation, using the explicit form of 
the interaction /, (99). We do not write down the analysis, one can follow 
closely e.g. [BFSS], [Ml], [BPS]. 

Lemma 5.1 We have the expansion 

PpinllPp = \Po{r + 0(e°))Po P(A < p) + O + , (182) 

where p is the parameter characterizing the infrared behaviour of the form 
factor (see Theorem 2.5; in the situation of Theorem 2.3 we set p = 0), 
0(e°) is an operator whose norm vanishes in the limit e — >• 0, and where 



[ da\g{l,a)\', (183) 
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and r is the non-negative operator on RanPo which has the following ma- 
trix representation in the basis {(po ® (fo, ■ ■ ■ , ^d-i ® ^d-i} of RanPo-' T is 
tridiagonal with diagonal given by [a, 1 + 2a, . . . , 1 + 2a, 1 + a] and constant 
subdiagonal and superdiagonal with entries —y^a{l + a), a = p(l) — -^Zi- 
The kernel ofT is spanned by the Gibbs state (52), ker(r) = Cf2i^^, and the 
spectrum of T has a gap j > at zero which is uniform in P > (3^, for 
fixed. 

Remark. It is easily verified by explicit calculation that fii^g is the unique 
element in the kernel of F. To see that the gap 7 is independent of (3 for 
large /3 one can use the fact that F converges to the matrix diag[0, 1, 1, . . . , 1] 
in the limit /3 — > 00. The latter matrix has a gap of size one. 
It follows from the lemma that 



29X^ ( PpIRliPp 



— + 4)). (184) 



where P^^p = l — Pi^^, and = is the projection onto the span 

of the Gibbs state (52). Using this estimate in (181) gives 

(5),, > min{^7,f}||^A||^-^7(A,/.P(A<p)),, 

Let us choose the parameters like this: e = A^^/^'^°, 9 = A^/^°°, p = X, 
p > —1/2. Then the minimum in (185) is given by (provided A < 

(47)-25/i3) and the error term in (185) is O (AVioo + o(x^)) = 0{X^). The 
virial theorem tells us that (-B)^^ = 0, so 

(Pi,^P(A<A))^^>l-0(A°). (186) 

We may write (186) as 

iPx = A,/3P(A < X)iJx + 0(A°) = Qi,f3 ® (P(A < X)xx) + 0(A°), (187) 

for some vector xx ^ T ®T with norm ||xa|| > 1 We point out that 

all estimtes are uniform in ^ in any compact set. This is easily seen by notic- 
ing that the only way ^ enters is through the term K^, which is uniformly 
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bounded in ^ belonging to any compact set in M^. This finishes the proof of 
Theorem 2.3. 

Proof of Corollary 2.4- We denote by Pi^/j, P/j^o and the projections 
onto the spans of ^i,/3, Jl^,o and il| see (52), (82) and (104). Since ||-P/3,o — 

x\\ ^ as X ^ (uniformly in in any compact set and in /3 > /3o, for 
any po fixed, [FM2]) it follows that 

P/3,oV'A + 0(A°) 

(Pi,/3®Pno)^A + PnoV'A + 0(A°) 
^^i,/3<8(Pqo^(A<A)xa)+0(A°) 

where we used (117) in the last step. It suffices now to observe that Pq^^P(A < 
A) converges strongly to zero, as A — > 0. This follows from = Pcijr <8) 
PQ^ + ll^®Pn^, 

P(A < A) = (P(dr(a;) < A) (g) P(dr(a;) < A))p(A < A) 

and the fact that dr(cj) has absolutely continuous spectrum covering 1R+ and 
a simple eigenvalue at zero, iljr being the eigenvector. ■ 
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